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I'^ACHIKE  SELECT JLK  STUDY 


INTRODUCTION 

The  selection  of  a  machine  from  a  class  of  available  machines  Is 
usually  a  question  of  some  Interest  and  complexity*  One  finds  available 
equipment  with  various  service  rates,  reliabilities  and  various  levels  of 
maintainability.  Assuming  that  purchase  price  or  rental  varies  across  the 
class  of  machines  it  is  usually  desired  to  select  the^^es^^Tiachine,  meas« 
ured  in  some  sense,  from  those  available  to  perform  a  specific  task*  -In  ■ 

*h3n  pT'**’  or,a  /sp  some  measures  of  performance  for  a 

very  restricted  class  of  machines*  One  might  expect  the  following  measures 
of  performance  to  be  important:  \ 

a*  Average  work  backlog. 

b*  Average  work  backlog  udien  machine  is  out  of  ccmmission* 

c*  In  commission  rate* 

d*  Machine  utilization  when  in  commission* 

e*  Fraction  'of  time  during  which  work  backlog  exceeds  a  specified  level* 
If  one  could  predict  these  measures  of  performance,  given  machine  relia¬ 
bility,  maintainability,  service  rate  and  arrival  rate  of  work  for  service, 

it  follows  that  the  effectiveness  of  a  pairticular  machine  might  be  assessed 

« 

and  choice  exercised*  To  illustrate  suppose  one  chose  average  woiic  backlog 
as  a  measure  of  machine  effectiveness*  Then  one  might  choose  that  machine 
which  generated  an  acceptable  bac||clog  at  least  overall  cost,  l*e*,  rental  atkd 
maintenance  costs*  Suppose  three  machines  were  available  to  perform  a  speci¬ 
fied  task  and  suppose  the  rate  of  flow  of  work  into  the  installation  were 
known.  One  might  find  that  the  following  chart  represented  the  average  back¬ 
log  versus  total  cost  for  a  fixed  arrival  rate  for  the  three  machines* 


<  % 


Fifi^re  1.01  -  Ilost  vs  Average  Jacklog 


Cost 


Therefore,  if  a  backlog  between  is  acceptable  then  machine 

3  is  the  best  choice.  If  a  larger  backlog  between  Bgjis 
acceptable  then  machine  2  is  best.  If  a  verj'^  large  backlog  is  accept¬ 
able  between  then  machine  1  is  best.  As  one  might  escpecty  no 

one  machine  of  the  three  is  best  for  the  whole  range  of  acceptable 
backlogs*  Further  figure  1  assumes  constant  average  arrival  rate  and 
would  undergo  alteration  for  a  'different  rate,  in  which  case  a  different 
machine  selection  might  be  involve'’  for  the  various  acceptable  backlogs* 

Clearlj'’  average  backlog  is  only  one  of  mar©'  ineasures  of  effective¬ 
ness  that  might  be  selects'^,  but  whatever  the  choice  it  is  usually 
possible  to  construct  a  chart  con;parable  to  figure  1  and  make  a  '•bestf* 
selection. 

The  paper  which  foUovrs  '^eals  only  with  the  characterisation  of 
candidate  equipments  in  terms  of  the  measures  of  performance  listed 
on  Page  1.  The  technique  fori-maklng  an  optimal  choice  will  appear 
in  a  subsequent  paper.  Farther,  the  paper  dpes  not  provide  a  ftilly 
developed  procedure  which  is  ready  for  implementation  and  use 
by  operating  personnel,  liather,  it  is  an  exploratory  study  whleh 
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»  K 


deala  with  a  h/pot^etloal  machine  and  process  >diLch  are  subject  to  severe  re- 
•trletlot>s.  The  studjr  is  presented  with  the  following  purposes  in  mind: 
a«  Tp  desNpmtrate  the  feasibility  of  the  approach. 

Ih  demonstrate  how^  although  technical  in  its  conception  and  develop- 
meat^  it  can  be  used  by  non-technical  personnel. 

c.  To  Indicate  the  course  vMch  has  to  be  taken  to  relax  the  restrlc- 

\ 

tlons  in  order  that  the  approach  may  become  more  general,  realistic  and  usable. 

The  paper  attempts  to  Introduce  the  non-matheroatlcal  reader  to  the  Ainda- 
mental  elements  of  probability  theory  ^Ich  fom  the  basic  structure  of  the  tech¬ 
nique.  It  is  believed  that  this  will  form  a  desirable  background  for  those  per¬ 
sonnel  idio  will  be  required  to  evaluate  the  approach  and  ultimately  decide 
whether  or  not  to  use  It. 

Fbr  those  udio  do  not  choose  to  follow  the  mathematical  rationale,  the 
final  portion  of  the  study  Includes  a  description  and  sample  of  the  nomographs 

and  tables  idilch  would  be  used  by  operating  personnel.  These  would  be  the  end- 

/ 

products  of  the  technique,  the  tools  tdiich  would  permit  operating  personnel  to 
predict  the  operation  of  the  various  candidate  equipments.  These  appear  on  pages 
21,  25,  and  29,  respectively.. 

1.  MODB. 

We  make  (perhaps  unrealistic)  restrictions  on  the  type  of  machine  and 
process  modeled,  in  the  Interest  of  indicating  a  type  of  approach  that  might 
be  taken  in  more  general  cases;  in  short  we  consider  this  an  eoqploratoxy  study. 

We  make  the  following  restrictions:  \ 

1.01  a.  Machine  capable  of  perfomlng  single  functions  (e.g.,  printer 

only  prints.) 

b.  ‘Machine  processes  only  one  unit  at  a  time  (no  parallel  processing). 

c.  Machine  functionally  disconnected  from  all  other  machines. 
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d.  Machine  service  time  much  less  than  one  hour  and  constant* 
e*  Machine  malfunctions  only  during  service  and  unit  In  service 

not  completed* 

f*  Arrival  of  units  requiring  service  Is  Independent  of  machine 

status* 

g*  Arrivals  non-lmpy  i*e*,  units  of  work  do  not  arrive  In  groups* 
h*  Order  of  service  is  first  come-first  served* 

We  may  summarize  by  saying  We  are  considering  machines  capable  of  one 
function*  rapid  In  performance*  that  process  units  sequentially  and  Independently 
of  other  machines,  Plurther  the  work  arrives  In  unite  rather  than  batches  and 
the  machine  processes  work  out  of  backlog  without  T>rlorlty*  The  machine  only 
malfunctions  during  operation  and  work  In  service  at  time  of  nwi  futxstlons  Is  not 
completed* 

2*  NA1URB  OF  TIME  DEUY  D15IRIHJT3DB3  CHOSEN 

We  assume  that  the  times  betWeen  malfunction^  between  repair^  and  be¬ 
tween  arrival  of  uidts  In  backlog  are  random  variables*  We  mean  by  this,  for 
example^  that  one  cannot  predict  exact  time  to  malfuhctlony  but  may  specify  the 
percentage  of  malfunctions  for  which  the  time  to  malfunction  will  exceed  a  par¬ 
ticular  value* 

For  example  in  Figure  1  we  have  a  particular  survivor  cuirve  for  a  machine* 
Figure  1  -  A  Survivor  Curve 


Fraction  of 
Malfunctions 

H(X) 


O  Xo  X| 

Time  to  Malfunction  (X) 
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The  survivor  curve  indicates  for  a  machine,  the  fraction  of  malfunctions  for 
which  the  time  to  malfunction  exceeds  a  specific  time.  This  curve  can  be  rep¬ 
resented  by  the  expression,  H(X),  where  X  expresses  time  in  some  chosen  unit. 

We  thus  see  that  the  fraction  of  malfunctions  for  vdiich  the  time  to 
malfunctions  exceeds  X^  is  .75,  Mathematically,  this  can  be  written  as 

2.01  PC  X^Xol=  H  (Xo)  = 

This  expression  states  that  the  probability  that  the  time  to  malfunction  (X)  is 
greater  than  or  equal  to  Xq  is  .75;  i.e,,  H(X),  at  the  point  Xq,  is  ,75.  Notice 
as  we  choose  a  larger  time  X^^  that  the  fraction  of  malfunctions  occurring  with 
time  exceeding  Xj^,  decreases.  Thus  if  Xj^^Xq  then 

2.02  H(Xi)*  PCX^xJ  4  PLX:?Xo3=H(Tfo). 

* 

When  X  =  0,  H(0)  =  1,  since  all  malfunctions  must  occur  subsequent  to  that  time. 
Now  suppose  that  one  observed  that  the  machine  had  not  malfunction  at  time  X^, 
Thus  the  machine  is  no  longer  a  so-called  (^zero  age^  machine  but  has  accumulated 
operating  time  Xq  since  its  last  malfunction  and  is  ^age^*  X^,  Since  we  know 
that  the  machine  must  malfunction  at  some  age  greater  than  or  equal  to  X^  it 
follows  that  H(Xq)  must  be  1,  (analogous  to  H(0)  as  described  earlier)  if  it 
is  to  properly  specify  the  survivor  distribution  for  an  '♦X^,  agedf*  machine.  If 
one  defines  a  new  function  >diich  gives  the  fraction  of  time  a  machine 

will  not  fall  in  the  Interval  of  time  X©  to  X^+Y,  given  that  it  has  survived 
to  age  Xq,  then  it  can  be  shown  that  the  new  function  will  be  related  to  the 
old  function,  H(X),  by  the  formula, 

u  iv'i- 

HCXo)  ^ 
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In  Figure  2.01  we  have  charted  H(X)  from  Figure  1.01  and 

Figure  2*01 


Ulus  from  2.03  we  see  that  IL  (Y)  Is  obtained  from  the  H(X}  curve  by  obtaining 

*0 

Its  value  at  Xq-^Y  and  dividing  by  H(Xq),  The  lmix)rtant  point  Is:  The  sturvlvor 
function  for  a  machine  of  any  age  may  be  computed  If  one  has  the  survlvDr  function 
for  a  ”zero  aged  machined. 

The  survivor  function  for  this  study  Is  assumed  to  have  the  form 

H(X)=  j  ©-70,  x:?o. 

We  have  in  Figure  2.02  a  graph  of  this  function 
Figure  2.02  -  H(‘X)s 


It  can  be  shoim  that  for  this  survivor  function  the  average  time  to  mal¬ 


function  is 


We  see  from  Figure  2*02  that  for  a  zero  aged  machine^  with  survivor  distribution 

-  '\ 

given  by  2,0U,  k5%  of  all  failures  will  occur  at  an  age  equal  to  or  exceeding' 
the  mean  tJjne  to  failure*  This  result  is  independent  of  the  value  for  the  mean 
failure  time. 

Now  suppose  an  item  were  known  to  have  age  Xq  and  had  expression  2,0k  for 
a  suivivor  function*  llion  we  see  that  the  survivor  functlcn  for  the  item^  given 
that  it  is  age  is 

a.o6  Hx^(V)= 

Notice  that  IL  (Y)  is  exactly  the  same  function  as  H(X)*  This  means  the 
sttrvlvor  function*  if  of  fom  2*04*  implies  that  the  accrual  of  age  on  a  device 
does  not  effect  its  chance  for  future  sxirvival*  In  other  words  the  likelihood 
that  an  equipment  will  fail  after  the  elapse  of  a  time  Y  is  not  related  to  the 
age  of  the  item*  It  turns  out  that  the  only  survivor  distribution  that  can  have 

t 

this  property  is  the  one  specified  in  2.0i!f* 

We  make  the  assumption  in  this  paper  that  malfunction  time,  repair  thne, 

« 

and  work  arrival  time  have  suivivor  distributions  of  the  form 


2*07 

HC-t)  = 

2*08 

H(>-)  = 

Bir)=  J  rZ’Ot&yo 

♦ 

2*09 

ri(a)  = 

;  eta)  -  f  j  ^ 

2.10 

E(t)  denotes 

• 

average  malfunction  time. 

2.11 

E(r)  denotes 

average  repair  time. 

-eiXo+W 
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2.12 


fi(a)  denotes  average  work  arrival  time, 


2.13  t  denotes  malfunction  time, 

2.14  r  denotes  repair  time, 

2.15  a  denotes  arrival  time. 

This  assumption  may  be  relaxed  if  actual  data  for  machines  indicate  a 
different  survivor  distribution  for  repair  and  malfunction  times,  or  if  items 
arrive  at  the  machines  with  a  different  survivor  function.  A  paper  given  by 
R.  M.  Lewis,  IBM,  at  a  conference  of  the  IRE,  1959«  titled  '*A  Study  on  the 
Collection  and  Analysis  of  Maintainability  Data,"  indicated  that  the  exponential 
distribution  characterised  repair  times  adequately  for  a  computing  system.  There¬ 
fore,  unless  evidence  to  the  contrary  can  be  offered  we  shall  maintain  the  above 
assumption  for  repair  time.  Exponential  distributions  have  also  been  used  ex¬ 
tensively  in  reliability  studies;  hence,  its  use  for  the  malfunction  distribution 

appears  reasonable.  While  the  exponential  distribution  may  be  unrealistic  for 

/ 

arrivals,  its  Impact  may  not  be  too  severe  for  heavy  workloads,  in  any  event, 
the  effect  of  relaxing  this  assumption  will  be  investigated  in  a  subsequent  paper. 

We  now  disciiss  the  relationship  between  E(t)  and  H(T).  To  do  this  we 
graph  in  Figiire  2.04  three  survivor  curves  for  the  three  different  mean  times. 
Suppose 

2.16  EaCt)  ^  . 

The  curves  labeled  1,  2,  and  3  are  associated  with  Ej^(t),  E^(t)  and  E^(t). 
Further  we  see  that  these  curves  intersect  in  only  one  point,  mainly  where  t  =  © 
and  H(0)  =  1. 
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Figure  2.04 


Ihree  Survivor  Curves 


Ihls  suggests  that  for  the  particular  curve  chosen  the  particular  H(t) 
is  uniquely  specified  by  E(t),  This  is  indeed  the  case.  Hotrever*  this  suggests  . 
another  fact;  that  if  one  knew  the  fraction  of  malfunctions  that  a  machine  sur¬ 
vived  for  a  specified  time,  T,  then  this  would  specify  the  Qiean  time  ^  failure, 
fi(t),  'IMs  is  in  fact  Irne  and  we  have 


ECt) 

T 


T  7  O  • 


T 

This  equation  permits  the  specification  of  E(t)  if  one  knows  for  a  specific 


time  T,  the  fraction  of  malfunctions  the  machine  accrues  for  age  larger  than  T, 


l.e.,  H(T).  Since  2.17  is  a  bit  Inconvenient  we  have  Figure  2.05  which  Is  a 

0 

graph  for  this  relationship. 
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Thus  if  H(t)  =  .7  and  if  T  is  5  hours  then  E(t)  is  13,89  hours.  In  other 
words  if  the  mean  time  to  failure  for  a  machine  is  13,89  hours  then  it  mil  oper¬ 
ate  without  malfunction  for  at  least  5  hours  for  70^  of  all  malfunctions.  Of 
coursoi  this  result  holds  only  if  the  survivor  curve  for  the  machine  is 

2,18  as  assumed  in  (2,07), 

Figure  2,05  can  be  used  to  estimate  E(t)  from  empirical  data,  Ih  order  to 
do  this  one  must  be  able  to  detemlne  a  value  for  H(T)  at  some  point  of  the  dis¬ 
tribution,  To  do  this,  the  following  procedure  can  be  followed: 

a.  Specify  a  time,  T, 

b.  For  a  large  number  of  malfunctions,  N,  count  the  number  that  occur 
when  the  machine  has  accimiulated  time  t  or  less  at  malfunction,  idiere  t  <  T,  Call 
this  number,  n, 

c.  Then  H(T)  =  N  -  n  ,  Using  Figure  2,05^  one  then  obtains  E(t), 

N 

or  Y,  Prom  this,  it  is  immediately  obvious  that  E(t)  “  YT, 

In  the  next  section  we  describe  the  operation  of  the  machine  and  process, 

t 

and  it  is  necessax7  to  speak  of  the  occurrence  of  a  machine  malfunction  in  one 
machine  service  interval  or,  in  short,  malfunction  during  one  service  time, 
FVirther,  we  often  speak  of  repair  completion  during  one  machine  service  Interval 
or  repair  completion  during  one  service  time.  Finally  we  refer  to  a  unit  work 
arrival  during  one  machine  service  Interval;  that  is,  an  arrival  during  one 
machine  service  time. 

Specifically  we  will  detezmlne  the  probability  that  a  machine  will  mal- 

I 

function  in  one  service  time  if  it  has  been  operating  for  time  Xg  since  the 
last  malfunction. 
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V/e  have  defined  in  2.03  the  survivor  distribution  for  a  "non-zero  age 


machine",  i.e.,  R.  (Y),  We  denote  a  machine  service  time  by 

*^0  r 


2.19 


X 


one  machine  service  time, 


^'to. 


Then  clearly  the  fraction  of  malfunctions  for  which  a  machine  aged  survives 

at  least  one  service  time  is  (A)»  Therefore,  the  probability  that  the 

o 

machine  does  malfunction  in  one  service  time,  given  that  it  has  been  operating 
for  time  since  the  last  malfunction,  is 

2.20  X|At,e=x»]= 


This  expression  states  that  the  probability  of  the  machine  malfunctioning  In  one 
service  time,  given  that  it  is  of  age,  X^,  at  the  start  of  the  service  period 
is  1  -  A).  Now  for  the  specific  survivor  function  used  in  this  study,  we 

have  from  2.06  and  2.07  that 


—  A 

2.21  Pc  x4  71 1 Ao]-  ^  5  }^70  . 


nils  result  states  the  very  important  fact  that  the  probability  that  the 
machinp  malfunctions  in  a  service  interval  is  independent  of  the  age  of  the 
machine. 

In  other  words,  if  the  survivor  distribution  is  as  specified  in  2.07  then 
the  chance  of  malfunction  in  a  service  interval  for  the  machine  is  not  affected 
by  the  operating  time  accumulated  since  last  repair.  It  can  be  ahown  that  if 
the  product"(^  is  small  then  (2.1^0)  is  approximately  equal  to  .  Written  „ 
mathematically,  we  have:  * 
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2.22  PC  0^  X4  ^  I  Koll“  ^  (vfhere  »  means 

approximately  equal).  The  error  committed  by  this  approximation  Is  less  than  or 
equal  to 


2.23 


Error  ^ 


Now  it  is  Important  to  note  that  the  mean  time  to  malfunction  was  specified 
in  2.07  as 


2.24  ECA)=  .  lliereforo  the  result  in  2.22  may  be  written 

in  an  interesting  form,  namely 


2.25  P  -  P[j>  4  y  ^ 


13^  =  i- 

Bix) 


ITius  the  probability  that  the  machine  will  malfunction  in  one  service  in¬ 
terval,  given  that  it  has  operated  for  a  time  is  apprcndmately  the  service 
time  divided  by  the  mean  time  to  failure  for  the  machine.  If  the  mean  failure 
time  for  the  machine  is  large  relative  to  the  service  time,  this  is  a  veiy  good 
approximation. 

Wo  lot 


2.26 

2.27 

2.26 


p.  = 

-  ^  A 

r-  e 

o<A  = 

7^ 

BU) 

rx  = 

1-  e 

Jl. 

eir) 

M  - 

1-  c 

- 

JZU 

Thus  and  ^  represent  the  probability  that  the  machine  will  mal- .. 

function  in  one  service  time;  the  probability  that  the  machine  repair  will  be 
completed  in  one  service  time;  the  probability  that  one  unit  of  work  will  arrive 
in  one  service  time,  respectively  and  these  are  only  a  function  of  average  failure 
time,  average  ropair  time,  average  inter-arrival  time  and  machine  service  time. 
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We  are  now  in  the  position  for 
3.  DESCRIPTION  OF  PROCESS' 

a.  Notation  employed, 

(1)  To  adequately  describe  the  process  by  which  the  model  was  deflYeid, 
it  is  necessary  to  define  several  of  the  basic  terms  used.  The  remaining  terms 
will  be  defined  later  in  the  discussion.  For  the  purpose  of  this  models  **i**  is 
used  to  indicate  the  status  of  the  machine  -  i.e,^  in  or  out  of  commission.  If 

i  =  li  the  machine  is  in  commission  and  conversely  if  i  =  0,  the  machine  is 
out  of  commission. 

(2)  is  used  to  indicate  the  magnitude  of  the  backlog  at  a  par¬ 
ticular  point  in  time, 

(3)  **t'*  is  used  to  indicate  time, 

b.  Verbal  description, 

(l)  Li  a  previous  section  of  this  report  we  have  specified  that 
average  rate  at  indiich  mrk  arrives  is  constant  over  time.  We  have  also  said 

that  not  more  than  one  unit  of  work  can  arrive  during  a  particular  service  time 

/ 

and  the  machine  can  only  service  one  unit  at  a  time,  Flirther^  the  machine  can 
go  out  of  commission  only  during  a  servicing  and  if  this  happens^  the  unit  in 
service  is  not  completed. 

Since  wo  have  made  the  above  restrictions  we  can  now  look  at  specific 
cases  of  backlog  growth  and  decay.  We  first  investigate  the  case  where  i  =  0 
(Machine  out  of  commission).  If  we  plot  the  backlog  J  against-  time  in  this  case 
we  get; 

I  Figure  3,01  -  Backlog,  liachine  out  of  Commission 


The  backlog  cannot  decrease  since  there  is  no  servicing  possible. 

We  now  investigate  the  case  when  i  •  1  and  the  backlog  is  initially  zero. 
Figure  3,02  -  Backlog,  Machine  in  Commission 


a  service  time,  the  backlog  with  the  machine  in  commission  cannot  exceed  1.  If 
the  backlog  does  exceed  1,  the  machine  must  then  be  in  the  state  previously 
described  when  i  =  0  (machine  goes  out  of  commission).  The  backlog  will  therefore 
grow  and  decay  fluctuating  between  0  and  1  as  long  as  the  machine  remains  in 
commission. 


(2)  Suppose  we  rK3w  investigate  all  the  possible  combinations  of 

t 

status  changes  as  time  advances  one  unit.  We  let  "t**  be  a  particular  instant 
in  time,  "t+l**  one  service  time  later,  ®  backlog  condition  with  the 

machine  in  commission,  a  backlog  condition  with  machine  out  of  commission. 

"a**  represents  an  arrival,  represents  no  arrival,  "jrf*  represents  a  malfunction, 
"jif*  represents  no  malfunction,  represents  a  repair,  and  finally  repre¬ 
sents  no  repair.  Therefore  the  following  table  can  be  constructed: 

Figure  3,03  -  Event  Tkble 


t/t+1 

CJ»oJ 

13-1,03 

CM.17 
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REH 

fJA] 

m 

a 

i 
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*  r 

0 

0 
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Let  us  now  look  at  a  few  examples  from  this  table.  Suppose  we  are  initially 
in  the  status  (J  items  in  backlog,  machine  in  commission).  To  move  to  a 

condition  in  the  next  interval  of  time,  we  would  have  had  to  have  no  arri¬ 
val  and  a  malfunction.  This  would  allow  the  backlog  to  remain  constant  but  the 
status  of  the  machine  would  go  from  1  to  0.  To  go  into  status  ^J-l,o]  would  be 
Impossible  (represented  by  0)  because  we  have  said  that  if  the  machine  malfunc¬ 
tions  (goes  from  1  to  0  status)  the  unit  being  serviced  will  not  be  completed 
and  therefore  the  backlog  can  not  decrease.  To  move  into  status  from 

,  we  need  both  an  arrival  and  a  malfunction.  Since  the  unit  being  ser¬ 
viced  is  not  completed,  the  arrival  increases  the  backlog  by  one  and  the  mal¬ 
function  changes  the  status  of  the  machine  from  1  to  0. 


Suppose  we  are  now  in  the  status 


[j,o]  (J 


items  in  backlog,  machine  out  of 


commission).  To  move  into  [M.SJ  or  [m,i]  is  impossible  because  with  the  ma¬ 
chine  out  of  commission,  no  servicing  is  possible  and  the  backlog  can  therefore 
not  decrease. 

The  remaining  cells  in  the  table  can  be  determined  in  the  same  manner  and 
the  completed  table  'will  thus  represent  all  possible  combinations  given  the 
restrictions  placed  on  the  model, 

(3)  Now  that  we  have  considered  all  possible  states  into  which  the 
machine  might  move  in  one  unit  of  tine,  we  may  graphically  represent  a  long 
range  status  sequence  for  the  pairticular  machine  in  question.  This  vrould  be 
done  in  much  the  same  manner  as  discussed  above  but  the  time  variable  would  be 
allowed  to  increase  over  a  long  period  of  time.  Hence: 


If  we  now  investigate  the  small  circled  section  of  this  graph  for  in¬ 
stance,  v/e  can  examine  in  more  detail  the  backlog  fluctriations  in  this  area. 
Figure  3.05  -  Machine  Backlog  (Short  Time  Span) 


In  this  chart,  each  interval  of  t  represents  one  service  time.  According  to 
the  restrictions  placed  on  the  process,  there  must  be  a  unit  in  backlog  at  the 
beginning  of  a  service  interval  in  order  for  a  service  to  take  place  during 
that  interval. 

If  we  look  at  the  Interval  from  t  =  0  to  1  we  see  that  the  backlog  was  0 
at  the  beginning  of  the  interval  so  there  is  no  service  performed  and  since  the 
backlog  is  still  zero  at  the  end  of  the  interval,  there  was  no  arrival.  The 
interval  from  1  to  2  also  has  no  service  but  the  backlog  Increases  from  0  to  1, 
hence  there  was  an* arrival. 

If  we  investigate  the  interval  from  5  to  6  we  see  that  the  backlog  has 
increased  by  one  indicating  that  the  machine  went  out  of  commission  (since  the 
backlog  at  the  beginning  of  the  interval >  O).  Looking  at  interval  7  to  8,  we 
see  that  such  a  condition  could  exist  in  one  of  two  ways  -  either  we  had  a  ser¬ 
vice  and  an  arrival  or  we  had  no  service  and  no  arrival.  If  we  consider  the 
fomer,  the  machine  goes  back  in^  comnisslon  at  point  7  and  if  we  consider  the 
latter,  at  point  8.  We  will  assimie  for  this  discussion  that  the  fomer  is  true. 

We  also  see  that  the  machine  again  goes  out  of  conmisslon  at  point  19  and 
comes  back  in  commission  at  either  21  or  22.  We  will  assume  it  comes  back  at  21. 

This  process  can  be  continued  for  an  indefinite  time  and  a  similar  pattern 
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to  the  above  could  be  expected. 

By  looking  at  the  long  term  process  (Figure  it  is  now  possible  to 

define  certain  terms  for  the  process, 

(4)  Definition  of  terms 

If  we  allow  the  intervals  in  which  the  machine  is  in 
commission  (0-5,  7-19,  etc.)  to  be  represented  by  t^^,  tg,  etc,,  to  t^  and  let 
the  total  time  elapsed  be  represented  by  T  we  can  define  the  fraction  of  time 
that  the  machine  is  in  commission,  P(l),  to  be 


3.01  p(i)=  - - 0.  =!^  ^ 

Conversely  by  letting  the  intervals  in  vrtiioh  the  machine  is  out  of  com¬ 
mission  be  represent^  by  a^  ag  ,,,,,, a^^  we  may  define  the  fraction  of  time  that 
the  machine  is  out  of  conmission  as  P(0)  so  that 

a.> 


3.02 


PCo)^ 


T 


Now  by  letting  ,,, 


measure  the  time  that  backlog  is  greater  than 


zei^)  in  each  interval  t^  through  t^^  (when  the  machine  is  in  commission)  we 
establish  a  definition  of  utilization  rate,  P[j>l/l3» 


3.03 


PLhilJ^ 


Ki  '■»  ♦'Kn 


S  k: 

ZiL 


(5)  Process  Averages 

Having  computed  the  foregoing  fractions,  it  is  now  possible 
to  compute  the  average  backlogs  for  the  process  >dien  the  machine  is  in  or  out  of 
commission,  and  we  are  also  able  to  compute  average  backlog,  over  both  conditions, 
which  one  might  call  the  unconditional  average  backlog. 
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If  we  let  bg,  •  •  •  average  backlog  in  intervale' vlflferi- 

the  machine  is  out  of  commission  we  can  compute  the  average  backlog  for  all 
Intervals  tdien  the  machine  is  out  of  commission^  Bp/oji  so  that 

0.1  bi-f  ba+*“-H2Hl»r*  __  ^ 


3.04 


EltJ/o]= 


Q.  I -f  4v+' + '  * '♦Q, 


Conversely  if  we  let  c^,  c^,  c^*  •  •  be  the  average  backlog  in  inter¬ 
vals  yihen  the  machine  is  in  commission  we  can  conpute  the  average  backlog  for 
all  Intervals  in  which  the  machine  is  in  coinni8slon«  t  ^bat 

Cltl  C.f.’tt. 


3.05  ECJI0  = 


■i|  *■  tv 


Xti 


Combining  the  two  above  formulae,  we  are  able  to  compute  the  unconditional 
backlog,  e|^j|,  so  that 

gdl(a£Q;-t-ED'IGX:t, 


3.06  ECJj* 


However 


a- 


3.07  Sai  v  =  T 

machine  is  either  in  or  out  of  coninlsslon,  not  both. 
Farther 


since  the 


pco;- 


PCI)  = 


T 


- 


by  3,01  and  3.02. 
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Therefore 


ELj]=  ECj|o1P(o)+  ED'lOPd) 


In  the  next  section  we  discuss  £|j/lJ,  P(0)  and  P(l)  as  a  function  of 
machine  and  process  characteristics^  1*6.,  the  average  failure  time^  average 
time  between  work  arrival,  average  repair  time  and  machine  service  time,  Ihese 
equations  were  obtained  by  establishing  a  system  of  difference  equations  which 
describe  mathematic all/  the  process  described  above  and  by  obtaining  the  gener¬ 
ating  function  for  the  process  from  these  difference  equations.  Once  we  obtain 
the  generating  function  the  various  averages  and  percentages  defined  above  may 
be  easily  obtained.  We  have  placed  the  system  of  difference  equations  and 
generating  function  for  the  process  In  the  appendixes  to  the  report, 

1,  RESULTS  OF  ANALYSIS,  In  Section  3  we  defined  severaj.  terms  for  the 
process.  In  this  section  we  specify  the  relation  between  these  terns  and  the 
parameters  employed  for  characterizing  the  machine  and  the  process  and  we  assume 
that  the  process  has  operated  for  a  long  enough  span  of  tjme  so  that  the  various 
measures  developed  have 'become  statistically  stable.  We  present  first  the  pre¬ 
cise  relationships,  then  simple  approximations  to  them,  provided  certain  assump¬ 
tions  are  made  and,  finally,  graphic  and  tabular  means  for  determining  the 
various  relationships, 

a.  Out  of  Corimisslon  Rate;  We  define  the  out  of  commission  rate 
to  be  the  fraction  of  total  time  that  the  machine  Is  In  reparable  status.  It 
can  be  shown  that  If  restrictions  1,01  are  In  force  for  the  machine  and  procesa 
that  the  out  of  conmlsslon  rate,  P(0),  Is 

^•01  R —  where  ^  and  P-  are  as  defined  by 

4  0-  fi) 

2,26,  2,27  and  2,28, 

f 
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A  simple  appraximatlon  I'or  this  result  may  be  obtained  If  the  mean  roalfctnsidon 
time,  mean  Interarrlval  time  and  mean  repair  times  are  large  In  relation  to 
the  machine  service  time.  It  can  be  shoim  that 


4.02 


Ed)*-  ,  It 


77  ^ 

E(f)» 

t(a.)7'7  ^  . 


than  a  service  Interval 


4.03 


Since  the  mean  malfunction  time  Is  much  larger 


_2_  _er2 

BtA.)  E(t)  ' 


We  therefore  see  that  If  mean  repair  time,  mean  malfunction  time  and  mean  Inter¬ 
arrlval  time  are  all  large  relative  to  a  machine  service  time,  that  the  fraction 
of  time  the  machine  will  be  out  of  commission  Is  approximately’  the  product  of  the 
ratio  of  service  tlhe  to  average  Interarrlval  time  and  mean  repair  time  to  mean 
malfunction  time.  If  this  is  ix>t  the  case  the  computation  for  P(0)  Is  facili¬ 
tated  by  employing  the  chart  in  Figure  4.01.  Flguire  4.01  Is  a  nomograph  which 
permits  the  determination  of  P(0)  with  virtually  no  computation.  It  would  be 
the  tool  actvially  used  in  the  implementation  of  the  technique  for  estimation  of 
the  In-oommlsslon  rate  of  a  candidate  equipnent. 

In  order  to  use  the  nomogra|>h,  one  starts  with  the  scale  In  the  center  of. 
the  lower  half  of  the  chart.  To  find  the  point  of  entry  on  the  scale,  take  the 
mean  time  between  arrivals,  expressed  In  hours,  £md  divide  It  by  the  service 
time  per  unit,  also  expressed  In  hours.  Find  the  point  on  the  scale  correspond¬ 
ing  to  this  ratio  and  draw  a  line  horizontally  to  the  right  until  it  reaches 
the  curve,  i.t  the  point  of  contact,  draw  a  vertical  line  until  it  intersects 


,  l'’.i.'ur'?  4.01  Cc-j.-pul^ition  cf  Ip-ocmmssicn  Rato , 

I  I  ■  I  I  I  I  I 
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one  of  the  lines  in  the  upper  right  hand  quadrant.  The  particular  line  to  be 
chosen  is  determined  by  the  ratio  of  the  mean  time  to  repairy  divided  by  the 
service  time^  each  expressed  in  hours,  and  further  divided  by  10,0CX},  If  the 
ratio  obtained  does  not  coincide  identically  vrith  one  of  the  lines,  visual 
interpolation  can  be  performed.  F^?om  the  point  of  intersection,  draw  a  line 
horizontally  to  the  left  xintil  it  intersects  one  of  the  lines  in  the  upper  left 
quadrant  of  the  chart.  This  line  is  determined  by  the  ratio  of  the  expected 
time  between  malfunctions  to  the  service  time,  in  hours,  and  further  divided 
by  10,000,  Again  interpolation  can  be  visually  achieved  if  the  ratio  obtained 
is  not  identical  to  the  value  shown  at  the  left  of  each  line.  At  the  point  of 
intersection,  draw  a  vertical  line  downward  tmtil  it  hits  the  horizontal  scale 
in  the  center  of  the  chart  to  the  left  of  the  ozdgin.  On  this  scale  can  be 
read  the  probability  that  the  machine  is  out  of  conmiselon.  Subtracting  this 
value  from  1  (one)  gives  the  ln>commission  rate. 

As  an  example,  assume  a  machine  has  the  follo^ng  characteristics: 

Service  .time:  200  lines  per  minute 

Workload:  Average  of  40,000  lines  per  S  hour  day 

Reliability:  Average  of  12  hours  between  malfunction 

Maintainability:  Average  of  7  hours  to  fix 

To  use  the  nomograph,  first  note  that  200  lines  per  minute  is  12,000  lines  per 

hour.  Average  service  time  per  unit  is  therefore  1/12000  hours  or  ,0000833 

hours.  A  workload  of  40000  lines  per  8  hour  day  is  an  average  of  5000  lines 

per  hour.  Average  time  between  arrivals  is  1  hours,  or  .0002  hours, 

.5000 

The  ratio  of  arrival  time  to  service  time,  i,e,,  E(a)  is  equal  to  2.4,  This 
is  the  figure  to  be  used  in  entering  the  chart  on  the  center,  lower  scale. 
Choice  of  the  proper  line  in  the  upper  right  quadrant  is  found  by  dividing  the 

% 


maintainability  figure  (7  hours)  i.y  the  service  time  as  Just  computed  (.0000833) 
and  again  by  10,000.  This  figure  is  8,4.  Ihe  appropriate  line  in  the  upper 
left  hand  quadrant  is  found  by  dividing  the  reliability  figure,  i.e,,  12  hours, 
by  the  service  time  (.0000833)  and  again  by  10,000,  i.e,,  14,4.  Use  of  the 
nomograph  then  shows  that  the  machine  will  be  out  of  commission  201"'  of  the  time. 
The  in-comrt.ission  rate  is  therefore,  80^, 

Ihe  machine  utilization  may  be  computed  by 


4.04  P[J''/I||]  = 


jy  M  Px 


where 


M  >  ji  j  fa  defined  as  before.  This  result  may  be  substantially 

simplified  if  the  mean  failure  ti/ne  and  mean  repair  time  are  large  relative  to 
the  service  time  for  the  machine. 

In  fact  we  have 


4.05 

*  mV 

"Jv  anc!  ECr)77/^,  Further  if  P(t)77F(r)  and  B(v)77^ 
then  we  have 


4.06 

4.07 


PCJXIO® 


I’  W 


e(r) 


M  . 


and 


Finally  IT  the  mean  interarrival  time,  3(a),  is  much  greater  than  the  service 
time,  ^  ,  we  have 


4.08 


PCdxIils. 

f 


2L- 

E(<!L) 
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Thus  we  see  that  simple  approximations  for  l/lj  are  quite  feasible 
if  the  mean  failure  time,  mean  repair  time  and  mean  interarrival  time  for  the 
process  are  all  large  relative  to  machine  service  time.  In  fact  if  mean  failure 
time  is  larger  than  mean  repair  time  and  if  mean  arrival  time  is  larger  than 
service  time  the  utilization  of  the  machine,  when  in  commission,  is  vezy  nearly 
the  ratio  of  service  time  to  mean  interarrival  time. 

On  the  other  hand,  rather  than  perform  these  computations  we  have  in 

I 

Figure  if«02  a  chart  which  pezniits  computation  of  PDwIO  graphically,  Ihls 
nomograph  is  the  second  tool  tddch  would  be  used  in  an  operational  environment 
and  would  pemlt  en^aluation  of  the  utilization  rate  with  virtually  no  computa¬ 
tion,  It  is  used  in  much  the  same  fashion  as  the  preceding  one.  The  center 
scale  on  the  lower  half  is  entered  at  the  same  point  as  in  the  first  nomograph. 

As  before,  a  horizontal  line  is  drawn  until  it  intersects  the  curve.  The  line  is 
then  extended  vertically  until  it  Intersects  the  appropriate  line  in  the  upper 

right  hand  quadrant  of  the  chart.  This  line  would  have  a  value  Identical  to 

/ 

the  line  used  in  the  left  hand  quadrant  of  the  first  nomograph.  Because  of  the 
nature  of  the  relationship,  however,  only  4  lines  are  shown  in  the  present  nomo¬ 
graph  with  values  of  1,  10,  100  and  over  100  respectively.  This  is  because  the 
computation  is  relatively  sensitive  to  values  of  B  (t)  between  1  and  10;  mur.h 

A 

less  sensitive  to  values  between  10  and  100;  and  virtually  insensitive  to  the  range 
of  numbers  beginning  at  100  and  increasing  without  limit.  As  before,  interpola¬ 
tion  between  the  lines  will  probably  be  necessary.  Prom  the  point  of  Intersection 

t 

(or  from  the  Interpolated  point  of  intersection),  a  horizontal  line  is  drawn  to 
the  left  until  it  intersects  the  appropriate  line  in  the  upper  left  hand  quadrant. 
This  line  is  determined  by  the  results  obtained  from  using  the  first  nomograph, 
namely  the  in-commlsslon  rate.  From  the  point  of  intersection,  a  vertical  line 
drawn  downwards  will  intersect  the  scale  at  the  utilization  rate. 
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Using  the  same  hypothetical  example  for  a  candidate  equj|}ment  as  we  did  In 
the  first  nomograph,  we  enter  the  scale  In  the  center  of  the  lower  half  of 
the  chart  at  2,k»  We  proceed  to  the  right  until  we  intersect  the  curve* 

At  the  point  of  intersection  we  move  upwards  until  we  hit  an  interpolated 
value  of  14*4*  We  then  move  left  to  the  line  with  value  of  *80,  the  In- 
commisslon  rate.  This  was  the  value  we  derived  from  the  first  nomograi^* 
Going  down  again,  we  read  k3%  as  the  utilization  rate. 

Thus  we  have  approximations  and  exact  functions  specifying  the  In- 
commlssjon  and  utilization  rates  as  a  function  of  mean  malfunction,  repair 
and  arrival  times  for  the  process.  We  turn  our  attention  to  the  average 
backlog  of  work  given  that  the  machine  is  out  of  commission,  l.e.,  E  [l/o]  . 


4.09 


E[j»o]=  I 


Ko)  l-M 

m  f.H-fi. 


1 


i^ere 


4.10  PoCO*  i*"  fj  "  ^ 

« 

A  slightly  more  convenient  form  is 


4.11 


E[J|oJ=  1  + 


Mil*  fife)- M*- 

feli-fi)  “Ml  fi+  fe) 


tdiere 


^  and  )i4.  are  defined  above. 

I  If  E(t)  and  E(r)  are  much  larger  than  ^  we  have 


^.12  EtJlol 


B(t)  E(r) 


M(t-U) 

^  A 


The  function  in  4.09  has  been  tabled  in  the  Efj/Oj  table  and  we  illustrate 
their  use  b,7  ';he  following. 
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Example;  Using  the  hypo'bhetlcal  piece  of  equipment  as  we  did  in  the 
nomographs^  we  calculate  and  /f  by  using  the  approximations  of 
2,260  2,27  and  2.28.  We  thus  have 


4.13  Rs 

=  .0000833 

12 

=  .00000694 

-  .694  X  10 

4.U  Pff 

-  .0000833 

=  .0000119 

s  .119  X  10 

E(r) 

7 

4.15  usr 

7 

=  ,0000833 

=  .417 

e^a)  ,0002 

The  table  is  divided  into  pages,  each  page  representing,  for  a  par¬ 


ticular  value  of  ,  the  average  backlogs  vdiich  would  be  expected  for  ranges 
of  ,  and  M,.  Thus,  one  first  finds  the  proper  page  by  matching  the 
compul>sd  as  in  4.14  with  the  at  the  top  of  each  page  of  the  table  and 
selecting  the  page  with  the  closest  ^  .  One  enters  the  table  by  choosing 
the  closest  and  ones  computed  as  in  4,13  and  4,15,  respectively. 

In  the  table  shown  as  Figure  4,03  each  entry  consists  of  two  numbers.  The 
first  number  is  .seven  and  the  second  two  digits.  We  notice  that  the  first 
seven  digit  number  has  a  decimal  point  to  the  left  of  the  first  digit,  Mstlce 
further  that  the  second  two  digit  number  has  either  a  minus  sign  or  no  sign 
following  it.  This  mode  of  representing  a  number  is  a  form  of  floating  point 
representation  and  is  read  as  follows. 

Rule;  (a)  If  the  sign  of  the  second  two  digit  number  is  positive,  move 
the  decimal  point  of  the  first  seven  digit  number  that  many 
places  to  the  right. 

Example;  ,2540073  04  means  2540,073 

(b)  If  the  sign  of  the  second  two  digit  number  is  negative,  move 
the  decirral  that  many  places  left. 

Example;  ,2500000  01-  means  ,02500000 
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We  notice  that  there  is  a  number  in  the  table  that  is  negative,  namely 
1010101-01.  This  number  does  not  represent  backlog  since  backlog  may'*t>ot 

be  negative.  The  nun.ber  means  that  if  one  attempts  operation  of  a  pro^ftfs 

/ 

with  the  parameters  indicated  bj'-  this  number  that  the  process  will  be  un¬ 
stable  and  no  steady  state  solution  exists.  Therefore  we  employ  1010101  -01 
to  indicate  process  parameters  idiich  create  instability  and  thus  little  or 
no  control.  If  one  attempts  operation  of  a  machine  at  the  parameter  levels 
Indicated,  and  if  the  assumptions  governing  the  study  hold,  chaos  will  result. 


We  notice  in  Figure  4.03  that  for 


II 

.9999999 

05-  ^  .0000199 

M-  = 

,4249998 

00  .4170000 

t|  - 

.5999999 

05-  .0000069 

that  E(j/0)  - 

.7636958 

05  i#e. 

E(V0)  = 

76369.58 

itens. 

This  figure  is  the  average  backlog  if  machine  is  out  of  commission.  Notice 

if  the  machine  service  rate  is  12,000  lines  per  hour  and  if  the  machine  goes 

back  into  commission  with  no  new  vrork  arriving  about  76370  =  6,4  hours 

12,000 

will  be  required,  on  the  average,  to  clean  up  the  backlog  of  work  on  hand. 

We  are  now  able  to  give  the  unconditional  average  backlog. 

It  can  be  shown  that  the  unconditional  average  backlog  is 

....  E.il-  . 


If  IS(r)  and  E(t)  are  much  larger  than  a  service  time  then 
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4.18 


CLV\  J 

EU)*  t-UlOJ  ^  E(t)+.ECv) 


Finally  the  average  backlog  when  machine  is  in  commlBslon  is 

We  have  developed  In-commlssion  and  utilization  rates  for  the  machine. 

Ihe  backlogs  for  the  process  under  various  machine  statuses  as  a  function 
of  average  malfunction,  repair  and  Interarrlval  time  and  machine  service 
time  have  also  been  developed.  Various  approximations  have  been  given  and 
nomographs  and  tables  developed  for  these  functions. 

5.  CONS  ms  ION.  It  is  thus  our  contention  that  a  machine  and  process 

may  be  characterized  mathematically  and  that  choice  among  alternative  machines 
may  be  quantitatively,  as  well  as  qualitatively,  exercised.  The  feasibility 
of  making  this  characterization  is  clearly  related  to  the  generality  one 
insists  must  be  incorporated  in  that  representation.  Hence,  we  have  shown 
that  if  one  restricts  a  process  severely  it  is  rather  easy  to  obtain  results 
which  permit  explicit  computation  of  process  measures  as  a  function  of  pro¬ 
cess  parameter.  The  question  of  validity  in  the  sense  that  such  a  characteri¬ 
zation  might  pemit  prediction  for  real  life  processes  is  another,  and  not 
altogether  unimportant,  matter.  It  would  seem  prudent  that  if  one  were  in¬ 
terested  in  "applicatiorf*  for  such  a  model  its  predictive  ability  should  be 
investigated  relative  to  a  real  life  process. 

Considering  for  a  moment  the  list  of  restrictions  enumerated  in  1.01 
it  seems,  on  the  basis  of  one*s  intuition,  that  the  most  severe  restriction 
in  the  list  is  "arrivals  be  non-lumpj'".  Experience  of  the  briefest  kind  in¬ 
dicates  othei  vise;  but  fortunatel.y,  this  restriction  is  easy  to  remove.  The 


restriction  to  a  single  machine  would  seem  equally  unrealistic  and  effort 
might  be  expended  in  relaxing  it*  It  is  believed  that  partial  relaxation 
of  this  assumption  can  also  be  obtained  without  too  much  difficulty*  The 
relaxation  of  these  two  assumptions  may  well  make  the  technique  workable  in 
the  sense  that  it  provides  reasonable  descriptions  of  real  world  activity* 

As  stated  above,  this  could  and  should  be  tested  against  real  data* 

should  be  relatively  clear  that  any  group  of  restrictions  in  the 
list  may  be  relaxed  but  usua3JLy  only  at  the  expense  of  complexity*  However, 
if  one  considers  the  benefits  of  extreme  generalization  worth  the  cost  then 
simulation  may  be  adopted  and  exceedingly  general  processes  studied* 

It  is  felt  that  consideration  of  machine  rental  and  service  rate,  to¬ 
gether  with  a  few  qualitative  observations,  are  not  adequate  in  characterizing 
the  process,  as  is  the  common  practice  at  the  present  time*  Indeed  the  machine 

rental  and  service  rate  assume  positions  of  importance  only  after  one  specifies 

$ 

the  acceptable  average  backlogs  and  even  then  the  mean  time  to  malfunction  and 

« 

mean  repair  time  are  perhaps  equally  important,  if  not  centrally  important. 

This  is  clearly  implied  by  the  backlog  equations  in  4.12,  A.l6  and  4*19* 
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